#市町村情報
shichoson_info <- st_drop_geometry(shichoson_geo)
gw_all_data <- st_transform(gw_all_data, st_crs(shichoson_geo))
intersections <- st_intersects(shichoson_geo, gw_all_data,sparse = FALSE) %>% which(. == TRUE)
shichoson_geo["passed"] <- 0
shichoson_geo[intersections,"passed"] <- 1
passed_pre_geo <- shichoson_geo[which(shichoson_geo$N03_001 %in% unique(shichoson_geo[intersections,1]$N03_001)),]
mapview(
passed_pre_geo,
zcol = "passed",
layer.name = NA,
col.regions = c("white", "red"),
popup = NA,
label = passed_pre_geo$N03_004,
legend = F
) +
mapview(gw_all_data,
label = NA)